Solar Physics 
DOI: 10.1007/' 



The Variation of Solar Wind Correlation Lengths 
Over Three Solar Cycles 



R. T. Wicks 1 • M. J. Owens 1 • T. S. 
Horbury 1 

© Springer •••• 

Abstract 

We present the results of a study of solar wind velocity and magnetic field 
correlation lengths over the last 35 years. The correlation length of the magnetic 
field magnitude A|g| increases on average by a factor of two at solar maxima 
compared to solar minima. The correlation lengths of the components of the 
magnetic field Xb xyz an d of the velocity \y YZ do not show this change and 
have similar values, indicating a continual turbulent correlation length of around 
1.4 x 10 6 km. We conclude that a linear relation between A151, VB 2 and Kp 
suggests that the former is related to the total magnetic energy in the solar wind 
and an estimate of the average size of geo-effective structures which is in turn 
proportional to VB 2 . By looking at the distribution of daily correlation lengths 
we show that the solar minimum values of X\gi correspond to the turbulent 
outer scale. A tail of larger Aibi values is present at solar maximum, causing the 
increase in mean value. 



1. Introduction 

The 22-year solar magnetic polarity cycle has been observed as the approxi- 
mately 11-year sunspot activity cycle for at least the last 250 years ( |Mursula, Usoskin, and Kovaltsov, 2001[ ).[ 
The effects of this cycle can be seen throughout the heliosphere: in the be- 
haviour of the solar wind (Richardso n and Kasper, 2008[ ), the flux of cosmic 



rays measured on Earth ( Forbush, 1954[ ), in auroral indices ( [Silverman, 19 92), 



geomagnetic activity ( |Cliver, B oriakoff, and Bounar, 19961, the nonlinearity of 



the response of the magnetosphere to the solar wind (Joh nson and Wing, 2005} , 



and in the rate of coronal mass ejections (CMEs) and hence the average \B\ at 



Earth (Owens and Crooker, 20061 



The interaction of solar particles and magnetic field with the rest of the 
heliosphere is mediated by the solar wind. The solar wind magnetic field and 



velocity fluctuations are turbulent (Goldstein, Roberts, and Matthaeus, 1995) 
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and the solar wind also contains many different forms of large scale struc- 
ture, such as CMEs ( |Wimmer-Schw cingru ber et al., 2006[ ) and co-rotating in- 
teraction regions (CIRs) formed by interacting streams with different veloci- 
ties ( Gosling and Pizz o, 1999 1. Turbulence and these larger structures affect the 
scattering of cosmic rays ( Minnie et al., 2003 1 and the interaction of the solar 
wind with planetary magnetospheres, causing auroral activity (Legrand and Simon, 1985). 

One method for studying turbulent fluids is to look at the turbulent correla- 
tion length A ( |Matthaeus et al., 1994[ ). This is the average distance over which 
correlation persists in the flow and is normally estimated from a single time series 
using the Taylor hypothesis ( |Taylor, 1938[ ). Direct measurements of the corre- 
lation length have been made using multiple spacecraft QMatthaeus et al., 2005 ; 



Wicks, Chapman, and Dendy, 20091, however they do not differ significantly 



from those made using single time series. Correlation length in the solar wind 
affects the scattering of energetic particles (Parhi et al., 20021 such as cosmic 
rays ( |Minnie et al., 2003 1 , it is an important parameter of models of the helio- 
sphere (fMatthac us et al., 1994| , giving an appropriate scale size for the models 
to resolve, and can be considered as a typical size for energy containing structures 
interacting with planetary magnetospheres. 

Given that the solar wind correlation length and the solar activity cycle 
both affect the scattering of cosmic rays and the behaviour of the Earth's 
magnetosphere it seems appropriate to study the solar cycle variation of the 
correlation length, its causes and its effects. In situ observations of the solar 
wind around 1 AU have been made regularly since the early 1970s using a 
variety of different spacecraft. Here we use the IMP8 and ISEE3 spacecraft, 
and the multi-spacecraft OMNI one minute resolution data set to provide mag- 
netic field and velocity observations from 1974 until 2009. This data set covers 
the last three solar cycles and so allows the current unusual solar minimum 



(McComas et al., 2008 Owens et al., 2008) to be compared to previous minima. 



We compare the correlation lengths to the changing sunspot number, CME rate, 
and magnetospheric activity as measured by the daily total of the Kp index over 
the most recent solar cycle, as well as looking at the general trends in correlation 
lengths for the last 35 years. 



2. Long term trends 

The data used solar wind velocity V and magnetic field B vectors and magni- 
tudes, all are one minute resolution. All observations are made close to 1 AU 
distance from the Sun. Here we use these data to provide a daily estimate of the 
correlation length of the velocity and magnetic field in the solar wind on time 
scales shorter than a day over the last three solar cycles. 

Correlation lengths are estimated by using up to 960 consecutive one minute 
cadence data points to calculate the autocorrelation function of the data, lagged 
by a time r between and 480 minutes, adding up to a total of 1440 minutes, 
allowing the calculation of a correlation length every day. An autocorrelation 
curve C(t) is calculated in this way every day and the time lag multiplied by 
the average solar wind speed is used to convert the time scale into a radial spatial 
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Figure 1. An example of the quarterly averaged spatial autocorrelation curves calculated 
for the solar wind magnetic field magnitude (solid line), GSE x component (dashed line), y 
component (dot-dashed line) and z component (dotted line). 

scale, that is C(x) = \V\C(t). The resulting spatial autocorrelation curve is a 
measure of the correlation of any scalar quantity (component of a vector) in 
the radial direction. This curve is then interpolated onto a standard spatial 
scale between and 9 x 10 6 km with a resolution of 3 x 10 4 km. This allows 
spatial correlation curves from different days to be compared directly. Finally 
a correlation length A is estimated from the correlation curves by fitting an 
exponential 

C(x)=exp(~). (1) 

Days arc only accepted if the spacecraft is in the solar wind and if there are 
over 480 out of 960 observations available. Individual days have very different 
correlation properties, caused by the very variable structure seen in the solar 
wind on timescales of a day such as CMEs, large amplitude Alfven waves and 
CIRs. 

To observe the change in A with solar cycle the daily autocorrelation functions 
are averaged over three months, giving a quarterly averaged correlation curve 
C(x) from which an estimate of A is made. An example of the resulting mean 
autocorrelation curves are shown in Figure Q] Quarterly averages are made only 
if there are at least 10 suitable days in the quarter of a year. The conditions on 
number of data points per day and number of good days per quarter are made 
to remove the statistically worst data points and keep a minimum standard 
between the correlation lengths calculated using different spacecraft despite the 
different orbits, data gaps and data qualities that occur. The IMP8 spacecraft 
had a geocentric orbit and so spent approximately half of its time in the solar 
wind. The OMNI data, however, are almost continuous, having very few missing 
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Figure 2. The quarterly averaged solar wind correlation length from 1974 until 2009. The 
top panel plots Ami, the second panel Xg XYZ . '^ ne third panel shows quarterly averaged 
sunspot number (R) as an indicator of solar activity. The fourth panel plots X\v\. The bottom 
panel shows the components of velocity Xvxyz > ^ v x 



is almost identical to A 



V|- 



days and fewer data gaps within individual days. Thus the quarterly correlation 
lengths since 1995 calculated from the OMNI data provide the most complete 
record of correlation over a solar cycle. The ISEE3 spacecraft data are used only 
during its primary mission (1979 - June 1982) before it was renamed ICE and 
left near-Earth space for comet Giacobini-Zinner. We have used the ISEE3 data 
specifically to fill in a period of time when the IMP8 data are particularly sparse. 

The quarterly correlation lengths are plotted in Figure [5J the top two panels 
show A calculated from the magnetic field magnitude \\gi and components 
\b X yzi the middle panel shows quarterly averaged sunspot number as an 
indicator of solar activity, and the bottom two panels show A calculated from 
the solar wind speed A|y| and velocity components Xv XYZ - The errors on the 
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estimations of A from these results are an order of magnitude smaller than the 
values themselves and hence are not plotted. Where the IMP8, ISEE3 and OMNI 
data sets overlap in time the correlation lengths calculated are similar in value. In 
particular between 1994 and 2001 both the IMP8 and OMNI data are available 
and although there is some systematic effect making the Omni correlation lengths 
slightly larger on average (particularly in the velocity measurements) they are 
not outside the scatter of the IMP8 results. This suggests that the estimates of 
the correlation length are consistent across all the data. 

We can also compare our values for correlation lengths to previ- 
ous estimates. Single spacecraft studies at 1 AU typically produce 
values of approximately Xb = 1.12 x 10 6 km and Ay = 2.83 x 10 6 km 
(Mattha eus and Goldstein, 1982 ). More recently it has become possible to 
use multiple spacecraft to calculate truly spatial correlations and the 
typical lengths of these are Xb = 1.19 x 10 6 km (Matthaeus et al., 2005). 
Wicks et al., (2009) made the first study of the change of multi 
point correlation length measurements with solar cycle and found 
that at solar minimum Ai_g| = 0.75 x 10 6 km and at solar maximum 
Ami = 1-75 x 10 6 km, but that the components of B had a roughly 
constant correlation length of Xb xyz ~ 1-1 x 10 6 km. These results, 
made using many different techniques, are broadly consistent with 
those shown here and make us confident that our results are reliable. 

The top panel of Figure shows that Ami changes with the solar cycle. The 
four activity minima contained in the data all give rise to similar correlation 
lengths of \\b\ ~ 2 x 10 6 km. For each of the last three solar activity maxima 
there is a peak in correlation length shortly after the peak in sunspot number, 
although the peak of solar cycle 21 is partially missing due to data gaps. The 
change in Ami between minimum and maximum is approximately a factor of 
two in all three cycles and the very highest values of A|_b| m each cycle are over 
three times the minimum values. The most recent estimates of X\b\ made here 
are from the end of 2008 and have not yet reached the previous minimum value 
and so it seems that solar minimum had not yet been reached by the end of 
2008, based on this measure. 

The second panel in Figure [5] plots the quarterly correlation lengths of the 
three GSE components of B. All three components have similar correlation 
lengths and remain approximately constant with average values of (Xb x ) = 
1.4 ±0.3 x 10 6 km, (X By ) = 1.4 ±0.3 x 10 6 km and (X Bz ) = 1.1 ±0.2 x 10 6 km 
over the whole 35 year time period. 

The quarterly correlation lengths of solar wind speed \V\ and the GSE V x 
component shown in the bottom two panels of Figure [2] are almost identical, as 
one would expect. The only clear trend in the behaviour of Aiyi is an increase 
in value since the most recent solar maximum. The previous highest values are 
all at the beginning of activity minima (1975, 1984, 1994), particularly around 
1975. The previous high values of Aiy-i from 1974 until 1980 could indicate that 
the current behaviour of the solar wind velocity is in fact a return to longer 
correlation lengths which previously occurred in the weaker cycle 20. However, 
without earlier observations than used here it is impossible to be certain. 



SOLA: R_Wicks_35_Year_Correlation_v4.tex; 11 February 2010; 11:52; p. 5 



Wicks et al. 



The correlation lengths of the GSE y and z components of V are plotted in the 
bottom panel of Figure [5J like the components of B they remain approximately 
constant for the whole 35 year period, with average values of (Xv Y ) = 1.7± 0.4 x 
10 6 km and (X Vz ) = 1.3 ± 0.5 x 10 6 km. The components of B and V (with the 
exception of V x since it is effectively \V\) have very similar correlation lengths 
and do not vary measurably with time. Thus the perpendicular components 
of V and all the components of B have similar correlation lengths. \\b\ at 
solar minimum is also approximately the same size. This is consistent with the 
presence of continual, developed Alfvcnic turbulence at 1 AU in the ecliptic, that 
is unaffected by solar activity. Further implications are discussed below. 



3. Solar Cycle 23 in detail 

Solar cycle 23 has the most complete data coverage of any of the three cycles 
studied here thus we investigate it further to identify the wider effects of the 
changing correlation length of \B\. To do this we compare the quarterly average 
\b\ with sunspot number, estimated CME rate, mean daily total magneto- 
spheric activity index Kp, and the quantity VB 2 over the last 14 years. These 
are plotted in Figure |3] similar plots have been made for the previous two cycles 
and the results are consistent, however solar cycle 23 has the best data coverage 
and so is the example used here. The data are split into minimum (•) and 
maximum (x) activity phases of the cycle, where maximum is defined as the 
quarters with an average sunspot number greater than 40. In order to 
compare correlation length with a consistent set of CME rates over the solar 
cycle, we take the same approach as ( |Owens et al., 2008 1, and combine the 



early duty-cycle corrected LAS CO CME rates QSt Cyr et al., 2000 ), processed 



LASCO CME catalogue rates (|Yashiro et al., 2004| | Gopalswamy et al, 2008] ) 



(See also http://cdaw.gsfc.nasa.gov/CME_list/index.html), and the recent duty- 



cycle corrected STEREO Corl CME rates made by St Cyr and Xie, available at 
|http: / /corl.gsfc.nasa.gov/| See ( |Owens et al., 2008] ) for more detail. 

The top left panel of Figure [3] plots the sunspot number against Ai^i. There is 
a weak trend showing increasing correlation length with sunspot number for the 
times defined as solar minimum, but the longest correlation lengths are found 
after sunspot maximum in the declining activity part of the cycle. Comparing 
A|b| with CME rate in the top right panel of Figure [3] gives very similar results 
to sunspot number, showing some relation but with considerable scatter. The 
bottom two panels of Figure [3] plot Ai b i against daily total magnetospheric activ- 
ity index Kp and the calculated solar wind quantity VB 2 , which is the amount 
of magnetic energy passing the spacecraft in a given time and a known driver 
for magnetospheric activity ( [Akasofu, 1979[ ) . Both show a linear agreement with 
the measured Aibi and reduced scatter compared to the CME rate and sunspot 
number. This supports the idea that A|_b| is related to the energy containing 
scale of the solar wind and that it is these \B\ scales that are geo-effective. 
Furthermore, this suggests that Kp or other magnetospheric activity indices 
could possibly be used as historical proxies for correlation in the solar wind 
magnetic field as well as its strength. 
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Figure 3. Plots of A|^| during cycle 23 against Sunspot number (top left), average daily CME 
rate (top right), magnctospheric activity index Kp (bottom left) and solar wind quantity VB 2 . 
The minimum phase of the solar cycle is plotted as • and the maximum phase as X to show 
the different behaviour during the different halves of the cycle. Maximum is denned as 
when the sunspot number is larger than 40. 



Finally we calculate the probability distribution of daily estimates of X\b\ 
during the two years with lowest average A|g| (1996 and 2008) and the year 
with the highest (2003) in the OMNI data set. These are plotted in Figure 
HI The daily correlation lengths are calculated by fitting the same exponential 
function to each day of data rather than to quarterly averages: the resulting 
correlation lengths are very variable but their distribution provides insight into 
what processes are causing the long term variation. The two minimum years 
have similar shaped curves with a peak at approximately 1.8 x 10 6 km and a 
slight tail towards longer correlation lengths. 2003 has a peak at a similar length 
to the minimum years but it is much lower and the tail of longer correlation 
lengths is much increased. 



4. Conclusions 

The averaged correlation lengths used here can be thought of as the average 
structure size in the variable being analysed. The values calculated here are con- 
sistent with previous estimates (Matthaeus et al., 2005; Wicks, ~Chapman, and Dendy, 2009] ) 
As such these results can be explained by the continual presence of turbulent 
structures with a constant average size of approximately 1.5 x 10 6 km throughout 
the solar cycle. This provides the minimum correlation length for \B\ and the 
peak in Figure 3] as well as the constant correlation length of the components of 
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Figure 4. Probability density functions of daily calculated correlation lengths from three 
years of the OMNI data set, two years with small correlation lengths close to solar 
minimum (1996, solid line and 2008, dot-dashed line) and the year with the largest av- 
erage correlation length, close to but not directly coinciding with solar maximum 

(2003, dashed line). 



B and V. This value is unaffected by the solar cycle and is consistent with the 
continual presence of developed Alfvenic turbulence at 1 AU. The correlation 
length of \B\ has an additional component added by solar activity causing it 
to change with the solar cycle, this could be partially caused by CMEs (Figure 
[3} but the highest correlation is found between Aibi and VB 2 , which is not 
necessarily related to any single type of event in the solar wind. 

The correlation of Ai^i with Kp and VB 2 suggests that the former is an 
estimate of the energy containing structure size in the solar wind. This could 
suggest that A|s| is related to the average size of packets of plasma in the 
solar wind which would provide a physical outer scale in slow wind in par- 
ticular. A|s| is also consistent with the idea described in QBorov sky, 2008) of 
large scale structures or flux tubes in \B\ which contain, or at least provide 
an upper scale for, the turbulent solar wind. This also implies that historical 
records of magnetospheric activity could be used to infer average structure size 
in the solar wind as well as the magnetic field strength and solar wind velocity 
dLockwood, Stamper, and Wild, 1999[ |Svalgaard and Oliver, 2007[ ). 

We note that A| v| has reached an unprecedentedly large value in recent years 
and that this might be connected with the current decrease in solar wind den- 
sity and ram pressure ( McComas et al., 2008[ ), although we currently have no 
explanation for this. X\b\ had not yet returned to its previous minimum value 
by the end of 2008 suggesting that solar cycle 23 had not reached minimum by 
this time. 
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